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TRANSVERSE VIBRATIONS OF HOLLOW THIN-WALLED 
CYLINDRICAL BEAMS 1 

By Bernard Budiansky and Edwin T. Kruszewski 


SUMMARY 

The variational principle, differential equations, and bound- 
ary conditions considered appropriate to the analysis of trans- 
verse vibrations of hollow thin-walled cylindrical beams are 
shown. General solutions for the modes and frequencies of 
cantilever and free-free cylindrical beams of arbitrary cross sec- 
tion but of uniform thickness are given. The combined influence 
of the secondary effects of transverse shear deformation, shear lag , 
and longitudinal inertia is shown in the form of curves for 
cylinders of rectangular cross section and uniform thickness. 
The contribution of each of the secondary effects to the total 
reduction in the actual frequency is also indicated. 

INTRODUCTION 

The elementary theory of bending vibration is often in- 
adequate for the accurate calculation of natural modes and 
frequencies of hollow, tliin-walled cylindrical beams. Such 
secondary effects as transverse shear deformation, shear lag, 
and longitudinal inertia, which are not considered in the 
elementary theory of lateral oscillations, can have appre- 
ciable influence, particularly on the higher modes and 
frequencies of vibration. The effects of transverse shear 
deformation and of rotary (rather than longitudinal) inertia 
have been studied by many on the basis of the original inves- 
tigations of Rayleigh (ref. 1) and Timoshenko (ref. 2). 
Anderson and Houbolt (ref. 3) have presented a procedure 
for including the effects of shear lag in the numerical calcu- 
lation of modes and frequencies of box beams of rectangular 
cross section. However, there does not appear to exist a 
general solution for the vibration of hollow beams that in- 
corporates the influence of all the secondary effects men- 
tioned. 

The purpose of the present report is threefold: First, to 
exhibit the variational principle, differential equations, and 
boundary conditions appropriate for the analysis of the 
uncoupled bending vibration of hollow thin-walled cylindrical 
beams; second, to give general solutions for cantilever and 
free-free cylinders of arbitrary cross section but of uniform 
thickness; and finally, to show quantitatively the influence 

1 Supersedes NACA TN 2682, “Transverse Vibrations of Hollow Thin-Walled Cylindrical 


of the secondary effects by means of numerical results for 
hollow beams of rectangular cross section of various lengths, 
widths, and depths. 


SYMBOLS 


A 

cross-sectional area 

A n 

Fourier coefficient 

■As 

effective shear-carrying area 

fit 

parameter defined in equation (30) 

r 

constant 

K 

modulus of elasticity 

G 

shear modulus of elasticity 

I 

moment of inertia 

K 

geometrical parameter defined in equation (29) 

L 

length of cantilever beam, half-length of free- 


free beam 

Ni 

parameter defined in equation (38) 

T 

maximum kinetic energy 

II 

maximum strain energy 

a 

half-depth of rectangular beam 

b 

half-width of rectangular beam 

Omni b n 

Fourier series coefficients 

i, j, m, n 

integers 

k B 

frequency coefficient, cu-y/ 

k s 

1 / E 1 

coefficient of shear rigidity, ^--y fl g 

km 

coefficient of rotary inertia, -^-y 

V 

perimeter of cross section 

s 

distance along periphery of cross section (see 


fig- 1) 

t 

wall thickness 

u(x,s ) 

longitudinal displacement in ^-direction 

w{x) 

vertical displacement in y-direction 

X 

longitudinal coordinate 

V 

vertical coordinate 

V 

y-coordinate of center of gravity of cross section 

7zs 

shear strain 
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e x longitudinal strain 

9 inclination of normal with vertical (see fig. 1) 

X Lagrangian multiplier 

p mass of beam per unit length 

p mass density of beam 

<r longitudinal direct stress 

r shear stress 

oi natural frequency of beam 

oiq natural frequency of beam calculated from ele- 

mentary beam theory 

5 t j Kronecker delta (1 if i~-j; 0 if i Aj) 

ip constraining relationship 

BASIC EQUATIONS 

Assumptions.— The problem to be considered is that of 
the natural bending vibration of a thin-walled hollow 
cylindrical beam whose cross section is symmetrical about at 
least one axis (see fig. 1). The transverse vibration is sup- 
posed to take place in the direction of this axis of symmetry 
of the cross section so that no torsional oscillations are 
induced. 

In the present analysis, the following simplifications are 
introduced: 

(a) Changes in the size and shape of the cross section are 
neglected. 

(b) Stress and strain are assumed to be uniform across the 
wall thickness. 

(c) The small effect of circumferential stress upon longi- 
tudinal strain is neglected. 

In accordance with statements (a) and (b), the distortions 
of the vibrating beam are completely described by the 
vertical displacement w(x ) of a cross section and the longi- 




tudinal displacement u(x,s) of each point of the median line 
of the beam wall. 

The longitudinal and shear strains are given in terms of 
u(x,s ) and w(x ) as 


bu 

tx ~bd 

(1) 

and 


bu , dw . . 

7zs ~bd + d^ Sm 9 

(2) 

and the corresponding stresses become 


TP 

° z==E bx 

(3) 

and 


n (bu . dw . A 
r x ,-6 U s + dx sin 6 ) 

(4) 


where 9 is the inclination of the normal with the vertical (see 
fig. 1). 

In elementary beam theory, where the effects of all shear 
distortion are neglected, the longitudinal distortion u(x,s ) is 
related to the vertical displacement w{x) by 


u (z,$)=(y-y) 


dw 

~dx 


where y is the y-coordinate of the center of gravity of the 
cross section. In the present report, however, u(x,s ) is 
allowed to be perfectly general, so that shear distortions (and 
consequently the so-called shear-lag and transverse-shear- 
deformation effects) are fully taken into account. Further- 
more, because cross sections are not constrained to remain 
plane, the inertia effect associated with motion in the 
longitudinal direction is more properly designated as the 
effect of longitudinal inertia than the effect of rotary inertia. 

Variational principle and geometrical boundary condi- 
tions. — The variational equation to be written is appropriate 
to beams whose ends are either fixed, simply supported, or 
free. For some such beam vibrating in a natural mode, the 
maximum strain energy is 

u -\5*f B Q£)' td,ix+ kKf G ( sf+S™ *)’*** 

(5) 



L 


X 


2 L 

I 

’ 

1 1 



(d) 1 

y 


(a) Typical cross section. (b) Sign conventions, 

(c) Cantilever beam. (d) Free-free beam. 

Figure 1. — Coordinate systems and sign conventions. 


where u{x,s ) and w{x) are the amplitudes of displacement for 
the particular mode considered. The maximum kinetic 
energy is 

!T=if (j> ptoi 2 w 2 ds dr+ij (|> ptoi 2 u 2 ds dx (6) 

where oi is the natural frequency of the mode under considera- 
tion and p is the mass density of the beam. The second 
term in equation (6) constitutes the contribution of longi- 
tudinal inertia to the kinetic energy. 

A natural mode of vibration must satisfy the variational 
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equation 

S(U~T)=0 (7) 

where the variation is taken independently with respect to 
u(x,s ) and w(x) and with the provision that both u(x,s) and 
w(x) must satisfy the geometrical boundary conditions of 
the problem; furthermore, u(x,s) must be periodic in the 
coordinate s with a period equal to the perimeter p. The 
geometrical boundary conditions are w=0 and u=0 at a 
fixed end and only w=0 at a simply supported end. At a 
free end no geometrical boundary conditions are imposed. 

Differential equations and natural boundary conditions. — 
Equations (5), (6), and (7) in conjunction with the usual 
procedure of the calculus of variations yield the following 
simultaneous integrodifferential equations for u and w: 


Fixed end: 

w=0 
u = 0 

Simply supported end: 

w=0 

El^= 0 
ox 

Free end: 

(j> Gi sin ej sin d ds — 0 

Et ^=0 
ox 

The integrodifferential equations (8) and (9), which specify 
equilibrium in the longitudinal and transverse directions 
respectively, can, of course, be written directly without 
recourse to the variational principle. 


(j > Gt sin sin 2 ds + M“ 2 w=0 


where 


n = (v pt ds 


( 9 ) 


( 10 ) 


and the boundary equations at each end of the beam are 


(10 




( 11 ) 


(^f Gt sin 0^ sin 6 ds Sw=0 (12) 


GENERAL SOLUTIONS FOR CYLINDERS OF UNIFORM WALL 

THICKNESS 

The following exact solutions for cylinders of uniform wall 
thickness are carried out by means of Fourier series in con- 
junction with the application of the variational condition 
(eq. (7)). This procedure, which does not require explicit 
consideration of the natural boundary conditions, was 
believed to be more expedient than a direct attack upon the 
simultaneous integrodifferential equations (8) and (9) and 
all their associated boundary conditions. 

Cantilever beam. — The geometrical boundary conditions, 
for a cantilever beam, as previously shown, are 


At a fixed end, both boundary equations (11) and (12) are 
satisfied by virtue of the fact that the geometrical boundary 
conditions require that both bu and Sw be zero. At a simply 
supported end 5ic=0, but, since Su(x,s) is perfectly arbitrary, 
the variational process forces the equality 

Et ^x=° (13) 

Finally, at a free end, since there are no geometrical con- 
straints, both bu and 8w are arbitrary and hence the varia- 
tional process forces, in addition to equation (13), the 
equality 

(j> Gt s * n s ’ n 0 ds=0 (14) 

Equations (13) and (14) constitute so-called “natural 
boundary conditions” because they are automatically satis- 
fied as the result of a variational process. Equation (13) is 
recognized as the condition of zero longitudinal direct stress 
while equation (14) simply stipulates that the total vertical 
shear force vanish. 

Thus to summarize, the appropriate boundary conditions 
required for the solution of equations (8) and (9) are 


ie(0) = u(0,s) = 0 

(see fig. 1). Appropriate assumptions for the displacements 
w(x) and u(x,s) are 

w(x) = C+ S b n cos~ (15) 

71=1,3,5 

and 

u(x,s)= S S «». Sin -fTj- cos — — (16) 

711=1,3,5 ?i=0,l,2 xXj p 

The condition «(0,s) = 0 is satisfied by each term of equa- 
tion (16); the condition 

w(0)=O+ 6„ = 0 (17) 

n = l,3 t 5 

is introduced into the variational procedure by means of the 
Lagrangian multiplier method. The choice of the particular 
trigonometric functions used in the Fourier series (15) and 
(16) was guided by consideration of the orthogonality 
required for the simplification of expressions in the strain 
energy. The constant C is needed in the expression for 
w(x) in order that w(L) be unrestricted. 
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Using equations (15) and (16) in equations (5) and (6) 
yields 

U-T= 

1 n,/ -A ^ mir m,TX 2nirs\ 2 j , 

O Et( X 2 »«« lyr cos Tf cos + 

2 Jo J \»»=1, 3, 5 71-0,1,2 V / 

ird>0*(ih s -««»- 

■i Jo J \m = l,3,5 n-0,1,2 


2 7i ir . niTrx . 2mrs 

sin — r sin 

p 2L p 


. ” , nir . mrx\- j , 

sine X; -Jj^rysin— s- ds rfj- 

91 — 1,3,5 ^-L// 


4 f Cf>«V f X & n 

-^Jo ./ \n = l,3,s 

I i; £ «« sin 

^ Jo J \ 7// = 1,3,5 71 = 0,1,2 


COS + ds dx- 


717 TTX 27l7rsV , , , to , 

— cos — ) ds ax (18) 
2L p / 


To make equation (18) stationary and at the same time 
satisfy the constraining relationship 


tp — C-\- X b„ — 0 

71 = 1,3,5 


it is sufficient to set 


S(U — T — X92) — 0 


(19) 

(20) 


where the variation is with respect to the a’s, b’s, and C 
considered as independent variables; here X is a Lagrangian 
multiplier. This variational process results in the following 
equations: 

b(U ~ U) _ a |£ = ± Gt * £ a u + 6 ) 2 ^ b f - 


dbi 


' dbi 


77 =0,1 ,2 


,2 L 2 


^-bt-i-lf? ~^ 2 C-\ 


lir 


= 0 


( 7 = 1 , 3 , 5 , . . .) ( 21 ) 

Gt^-Ajbi-Jpt^i 1 +8 w )o„=0 

(?'= 1 , 3 , 5 , . . .) 

(.7 = 0, 1,2, . . .) (22) 

D(U-T) .d<p 2 “ 2 L ,X'; no x 

"He X dC~~^ n£i«7 ( ) " " 1X0310 


where 


= 0 
A 


Hf* 

-f 


. . . 2nirx , 

sin 9 sin as 

V 


t sin 2 8 ds 


(23) 

(24) 

(25) 


With the use of the nondimensional parameters 


‘•■=17 

(26) 

, 2 El 
Cs A S GL 2 

(27) 

, , / I 

Cm ~ptL*~AL* 

(28) 

16 1 

Astf 

(29) 

B 2 =i 2 -1c R , 2 k n 2 (- ) 

(30) 


and 


equations (21), (22), and (23) may be reduced to 


“ inir 2 Lt 1 

X ~Ab 2 A r> 

71=0,1,2 4ft: S As 6 


(iir\ 2 1 

V27 ks 


LA 2 , 2/7 i 3 X . 

(- 1 ) 2 -^ kB °-]n =0 


k Ii 2 b i 


(i= 1,3,5,...) (31) 


(k s 2 Bi 2 +K 2 j 2 ) (1 -Ho,) a {J +K 2 Ajijbt—0 

(i=l,3, 5, .. .) 


k B 2 X — (-i)V6»+^ 2 6’+~J v = 0 

71 = 1, 3, 5 Tlx ' till 


(3= 0,1,2,...) (32) 
(33) 


For j=0, equation (32) becomes 


k.J i 2 -k,/k [{1 2 0 J a i0 = 


= 0 (»= 1,3, 5, . . . ) (34) 


Equation (34) is not coupled to any of equations (31) to 
(33). A given value of a i0 corresponds to the amplitude of 
the fth mode of longitudinal oscillation, and if this value of 
a i0 is not equal to 0, then equation (34) simply gives the 
frequency of this longitudinal mode. Consequently those 
equations in equation (32) for values of j—0 are not asso- 
ciated with transverse bending and so are ignored henceforth. 
For the remaining values of j (that is, ,?H0) equation (32) 
yields 


-K 2 l L Ajij 

a,i= ~k;w+K 2 7 bt 


(7=1, 3, 5, . . .) 

(j= 1,2,3,...) (35) 


Substituting the expression for a i} in equation (35) into 
equation (31) and solving for b ( gives 


b t - 


/ ,^2 , I 3 X 

N. 


(7=1,3, 5, . . . ) (36) 
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where 


N<=- 


v i 

8 k s > / j 


2_ 7 r 2 K 2 -4 . 2 ,. , 2 

16 1 






ra=l, 2, 3 


In the appendix this expression for N t is shown to be equivalent to 


Ni 32 B * 16 A s Bi ' ks2 Zj KViks^B^+K 2 ^) 2 1 


iV A 

16 A s 


2 

n = 1,2,3 


AJ 


(37) 


(38) 


Since the series in equation (38) is considerably more quickly convergent than that in equation (37), equation (38) should be 
used in actual numerical calculations of N t . *r 

Substitution of equation (36) into equation (33) and the constraining-relationship equation (19) gives the following two 
homogeneous equations in C and X: 


fi+£ s 2 ± £ (- 1 )” 2 

L 71 = 1,3,5 \7iTT / IXnJ L 71 = 1,3,5 


_2 1_ - 

^ = 0 

(39a) 

n-ir N„_ 

El 

^=0 
El u 


(39b) 


Finally the condition for a nontrivial solution for C and X gives the frequency equation 




03 n — l n -I 

H-W s (-1) 2 

n= 1,3,5 TlTT IV n 


«* 9 l 

1+^ B 2 s (-1) 2 fzir 

W = l,3,5 71 TT IX n 


Zi W 

n = l,3,5 iVn 


= 0 


(40) 


which the frequency parameter k B must satisfy. Since the 
terms of the infinite series which appear in the frequency 
equation contain k B itself, the roots of equation (40) are most 
conveniently found by trial. Fortunately the infinite series 
in equation (40) as well as the series in the definition of A T t 
converge rapidly so that only a few terms are needed to 
evaluate them with sufficient accuracy. 

Once k B has been determined for a particular mode, the 
corresponding mode shape can be found by letting C= 1 and 
solving either of equations (39) for X and then finally evaluat- 
ing b t and a tJ successively from equations (36) and (35). 

Free-free beam — symmetrical modes. — If the origin of a 
free-free beam of length 2 L is taken at the midspan (see fig. 
1), the form of the Fourier series assumed for w(x) and u(x,s) 
when the beam is undergoing a symmetrical mode of vibra- 
tion may be exactly the same as that assumed for the 
cantilever beam of length L (see eqs. (15) and (16)). The 
only difference in the ensuing calculations is that the con- 
straining condition (19) is not introduced. Consequently, it 
can be readily seen that the frequency equation for the 


symmetrically vibrating free-free beam is obtained from 
equation (39a) by setting X = 0 and is 

k B 2 [i+k B 2 £ (Y) 2 1= 0 (41) 

[_ 71 = 1,3,5 \7lwJ A n ._| 

After a particular root k B is found from equation (41), the 
shape of the corresponding symmetrical free-free mode may 
be obtained from equations (36) (with X = 0) and equations 
(35). 

Free-free beam — antisymmetrical modes. — Consider a 
free-free beam of length 2 L undergoing antisymmetrical 
vibrations. Explicit consideration need be given only to the 
right half of the beam (see fig. 1), and for this half-beam the 
only geometrical boundary condition that must be imposed 
is that w(0) = 0. The spanwise displacement u{ 0,s) is unre- 
strained by virtue of antisymmetry. 

Appropriate assumptions for the displacements w(x ) and 
u(x,s) are then 

w(x)= £ b n sm^+Cx (42) 

n=2,4,6 
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» » m-icx 2mrs , 

U{X,S)= 2j X a mn COS COS (43) 

m*= 0,2,4 n= 1,2,3 P 

The linear portion C'r of the expression for w(x) is needed in order to give the beam sufficient freedom at the tip (x=L). 
The choice of the particular trigonometric function in the series expansion for u(x,s ) was, as in the case of the cantilever 
beam, guided by consideration of the orthogonality required for the simplification of the expressions in the strain energy. 
The zeroeth term in the series for u(x,s) in the s-direction was omitted because it only leads to the frequency equation for 
longitudinal oscillations. 

Using equations (42) and (43) in equations (5) and (6) yields 


■T-=U' L <£ Et( ± X sin Tr c° s -T*y ds d-x+ 1 9 p(f> Gt 2 ± - 

6. In J \m=0,2,4 n = l,2,3 2X, 2A p / ZJ 0 J L» = 0,2,4 » = 1,2,3 

sin tf^X 0 K cos ^+CJJ ds dx-±£(£ P t<A * b n sin ^+rxj ds d 


2mr mir.r, . 2nirs , 

l,nn cos or Sin \- 

p 2L p 


\ pt " 2 ( X X « 

4Jo J \m = 0,2,4 n = 1,2,3 


2nirs\ 2 


mn cos 2 Z ™s — 


The variation of equation (44) with respect to the a’s, I where is defined in equation (37). 


b’s, and C gives, after suitable simplification, 

(B ( 2 k s *+K*j*)a (j -K 2 * Ajijb j=0 (i= 2,4,6, . . .) 

(3= 1,2,3, - - .) (45) 


Substitution of equations (50) and (51) into equation (48) 
and simplification gives as the frequency equation for the 
antisymmetrically vibrating free-free beam 


j 2 — k s 2 k UI 2 k B 2 


a 0j -K 2 ^ A i jC= 0 

0= 1,2,3, . . .) (46) 


k k dL+l A 

B L^4,6 \nir) N^* 2 As 


km 2 X 


B 0 2 k s 2 +K 2 n 2 


2 ^/( 2 ) 6i + 2 fcs26i 


(-l) i/2 -f- k B 2 CL=0 (i = 2,4,6,...) (47) 

XlT 


i 1 It . uz 

1 , 2 a n7r -A- n a 0n > 2 

n= 1,2,3 nts Vg 


(^'= 1 , 2 , 3 , . . .) (46) After a particular value of Ar B is found from equation (52), 

the shape of the corresponding antisymmetrical free-free 
1 If i« 7 r 2 . 11 / iir\ 2 , , _1 , 2 l _ mode may be obtained by giving C the arbitrary value of 

m^ 2,3 k s 2 A s 4 " a ' n 2 k s 2 V 2 / 1 2 B 1 unity and calculating the b’s and a’s successively from 

equations (51), (50), and (49). 

iy /2 _2_ kg 2 CL=0 (i = 2 4 6 , ) ( 47 ) Discussion of parameters. — The parameters entering in 

l7r the frequency equations merit discussion. The unknown 

^ £ QL natural frequency is contained only in the frequency co- 

T~2 ~T n ' 1T A n a 0n — ’ ~J—~2 — / J£T 

11 = 1 , 2,3 r4 s fc s efficient k B , which is defined bv the formula co~k B -\ — ^ 

V 

^ 2 jp 2l t—\yn ft _|_ A k ?CL — 0 (48) aIK ^ * s * n c °minon use in beam-vibration analysis. The 

8 Kx n 3 B parameters Ic s and k Rr are identical with the shear and 

inertia parameters defined in reference 4, which considers 
From equation (45) the effect of only transverse shear and rotary inertia on 

2 p beam vibrations. The quantity A s which appears in the 

Ajij present definition of k s is actually the effective shear-carrying 

a ’ i= B 2 k. A A K 2 ] 1 b ‘ ('f — 2,4, 6 , • . . ) (49) area when plane sections are constrained to remain plane; 

(. 7 = 1 , 2 , 3, . . .) that is, when shear lag is neglected. The remaining param- 

, . „ , . . eters appearing in the present derivation, namely, A/As, 

which, except tor sign, is the same expression as that ob- R and ^ ^ are essentially shape parameters which 

tamed for the cantilever and symmetrically vibrating free- actuaUy depend only 

on the contour of the cross section; 


k B 2 X ^~{-l) nl2 b n + \k B 2 CL = 0 (48) 

71 = 2, 4,(5 72-7T O 


From equation (45) 


K 2 {y Ajij 


l iJ TJ 21 2_|_74 r 2-2 (f 2,4,6, • ■ . ) (49) 

B t k & +Kj (j=l, 2, 3, ... ) 


tamed for the cantilever and symmetrically vibrating free- actuaUy depend 0 nly on the 


free beams (eq. (35)). From equation (46) 


R2 1 A « 

a oj B 0 2 k s 2 +K 2 j 2 ^ ^ 1,2,3, . . .) (50) 

Substitution of equation (49) into equation (47) gives 

bi= ~ { ~ l)m tw i CL ^= 2 ’ 4 ’ 6 ’ ■•.) ( 51 ) 


as shown in the appendix, 


As jl , 2 

A O L ' CXn 


A 2 re fif 2 ,3 


T7-, 2 A A „ 2 

T 2 A Sn ^ n 2 


and the /l„’s are simply the Fourier coefficients of the 
function sin d, which is dependent only on the shape of the 
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cross section. These shape parameters are related to shear- 
lag effects and their interaction with transverse shear 
and longitudinal inertia. 

The effect of longitudinal inertia is associated with the 
parameter km- If the effect of longitudinal inertia is to be 
neglected, it is sufficient to set k RI equal to zero in the final 
frequency equation. If k B! is equal to zero, B t becomes in- 
dependent of k B . Appreciable simplification in a trial-and- 
error solution for the natural frequency then results since, 
with B t independent of k B , the infinite summation contained 
in Ni is also independent of k B and need be calculated only 
once for any particular beam. As is shown in the following 
section, the effect of disregarding the influence of longitudinal 
inertia may often be negligible. 

Without presentation of details, it may be mentioned that 
for the case of a circular cylinder, which has no shear lag, 
all the /l„’s except Ai vanish and the frequency equations 
(40), (41), and (52) may be put into closed forms identical 
to those given in reference 4. Again, if in the general fre- 
quency equations k s is set equal to zero, the equations may 
be put into closed forms equivalent to those of reference 4 
where only rotary inertia is considered. 

RESULTS FOR CYLINDRICAL BEAMS OF RECTANGULAR 
CROSS SECTION 

In order to show quantitatively the effects of shear lag, 
transverse shear deformation, and longitudinal inertia on 
the natural frequencies of hollow tliin-walled cylindrical 


beams, numerical calculations have been performed for 
cylinders of rectangular cross section oscillating as free-free 
beams. The calculations have been limited to sy mm etrical 
modes of vibration, and consequently the frequency equation 
(41) is applicable. For rectangular cross sections the quan- 
tity Nf mat 7 be put into closed form as shown in the appendix, 
and this closed-form version of N t was used in the calcula- 
tions. A value of E/G equal to 2.65 (appropriate for 
aluminum alloys) was assumed. 

The results of these calculations are shown in figures 2, 3, 
and 4. In figure 2, the ratio of the natural frequency o> to 
the natural frequency co„ obtained from elementary beam 
theory is shown as a function of the plan-form aspect ratio 
L/b for cross-sectional aspect ratios of 1.0, 3.6, and <». 
The contribution of each of the secondary effects to the total 
reduction in the natural frequency for the cross-sectional 

aspect ratios ^=3.6 and 1.0 can be seen in figures 3 and 4, 

respectively. The cross-sectional aspect ratio of — = °° 

corresponds to the limiting case of a beam where the effects 
of transverse shear deformation and longitudinal inertia are 
negligible and therefore the reduction in natural frequency 
is due entirely to shear lag. 

The dashed lines in figures 3 and 4 show the reduction in 
frequency due to the inclusion of the effect of only transverse 
shear deformation as obtained from reference 4. 

The long- and short-dash lines are calculated from the 
frequency equation (41) with k RI = 0 and consequently 



(a) First symmetrical mode. (b) Second symmetrical mode. (c) Third symmetrical mode. 

Figure 2. — Change in the natural frequency of a symmetrically vibrating free-free cylinder due to the inclusion of secondary effects. 
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L/b 


Figure 3.— Contribution of transverse shear deformation, shear lag, 
and longitudinal inertia to the reduction in natural frequency for 



1.0 
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.6 


.2 



10 


14 
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Figure 4. — Contribution of transverse shear deformation, shear lag, 
and longitudinal inertia to the reduction in natural frequency for 

-=1.0. 

a 


represent the reduction in natural frequency when both shear 
lag and transverse shear deformation are taken into account. 
Thus the hatched area between the dashed and the long- and 
short-dash lines may be considered as showing the additional 
reduction in natural frequency when the influence of shear 
lag is considered. Finally, the solid lines are calculated with 
lc RI taken into account, and consequently the shaded area 
shows the additional influence of longitudinal inertia in 
reducing the frequency. 

Examination of figures 3 and 4 and the curves for — = °° in 

Cb 

figure 2 shows that the influence of shear lag increases as the 


cross-sectional aspect ratio increases; whereas the influence 
of transverse shear and longitudinal inertia decreases with 
increasing cross-sectional aspect ratio. Indeed, it appears 

from the results for — =3.6 that for this aspect ratio the 
a 

effects of longitudinal inertia may already be considered 
practically negligible. 

A word of caution concerning the interpretation of figures 
3 and 4 may be in order. Since in some cases the depth of 
the hatching increases with increasing L/b, it might appear, at 
first glance, that the shear-lag effect increases with increasing 
plan-form aspect ratio. However, if the additional effects 
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of shear lag are considered on a percentage basis with the 
dashed line as a base, it will be found that shear-lag effects 
actually reduce in percentage with increasing i/6. A s im ilar 
criterion should be used in judging the influence of longi- 
tudinal inertia. 

CONCLUDING REMARKS 

The numerical calculations show that secondary effects 
have appreciable influence on the natural frequencies of 
rectangular box beams of uniform wall thickness. These 
results constitute an indication of the probable inadequacy of 
elementary beam theory for the vibration analysis of actual 
aircraft structures of the monocoque and semimonocoque 


type and emphasize the need for practical calculation pro- 
cedures for such structures that would take into account 
transverse shear deformation, shear lag, and, when necessary, 
longitudinal inertia. The general solutions presented for 
cylinders of uniform thickness, as well as the numerical results 
for rectangular box beams, should be useful in the assessment 
of the accuracy of any procedure of this kind that may be 
developed. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Ya., January 21, 1952. 


APPENDIX 


Expressions for A s /A, I, and K 2 . — If sin 6 is expanded into 
a Fourier series 

2nirs 


TRANSFORMATION OF PARAMETERS 

or 


sin 9— Y) d, sin 

71 = 1 X 3 V 


(Al) 


the Fourier coefficients A n are the same as those defined in 
equation (24) ; that is, 




■ . 2nirs , 

sin 9 sm as 

V 


(A 2) 


The effective shear area A s (eq. (25)) can now be written 
as a function of the Fourier series expansion for sin 9 as 

A a ^ j > t [s A Sin ~^Jds (A3) 

With the use of the appropriate orthogonality conditions, 
equation (A3) becomes, after the integration is performed, 

s A n 2 =4 S A* 2 


« = 1 ( 2,3 


77 = 1,2,3 


or 


-As ,J_ nS a 2 

n 2—1 


A 2 71 = 1,2,3 


(A 4) 


The moment of inertia / of a cylinder is defined as (see 

fig- 1) 

/= yH ds — Ay 2 (A 5) 

where y is the y-distance to the center of gravity of the cross 
section and is given by 


But 


( P ytds 

y= JjL ~r- 

pt 

j— I sin 9ds 

Jo 


(A6) 
(A 7) 


V 2nir V 


2mrs \ 
P ) 


and, consequently, 


V= S A n A2- 

71 = 1,2,3 ZRir 


(AS) 


(A9) 


With the use of equations (A8) and (A9), the expression 
for I in equation (A5) becomes 


Ap 2 


AJ 


Q 2 ^—1 „ 2 

OF 71 = 1,2,3 n 


(A10) 


With the series expansion for I in equation (AlO), the param- 
eter K 2 , as defined in equation (29), becomes 


o A co a 2 

k 2 = 4 #- s — r 

ir A s 71=1)2,3 n* 


(All) 


Transformation of expression for N t . — In equation (37) 


N t was defined as 


N t 


iV A 


16/c.v 


s 


n 2 A„ 2 


As^ika'Bf+KW 2 


o^'a 2 (A 12) 


The infinite series that appears in this expression converges 
as A n 2 and therefore is a relatively slowly converging series. 
In order to increase its rate of convergence, the following 
transformations are made. 

By adding and subtracting A n 2 /K 2 inside the infinite 
summation in equation (Al2) and using equation (A4), the 
equation simplifies to 

N <-W b 'M, L-aAW -5 *•’ < A13 > 

By adding and subtracting A n 2 /K 2 n 2 inside the infinite 
summation in equation (Al3) and using equation (All), the 
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expression for N t can be transformed to 


Nt ~ 32 B ' 16 A a Bi ' la \ 1 ^tEhiKk*'Bt+KW) 2 k » (Al4) 

The infinite series in equation (Al4) converges as AJ/n* and therefore is considerably more quickly convergent than the 
series in equations (A12) and (Al3), which converge as A n 2 and A n 2 /n 2 , respectively. 

Closed form of N { for cylindrical beams of rectangular cross section, — For a cylindrical beam of rectangular cross 
section, with dimensions as shown in figure 2, it is possible to write the expression for N t in a closed form. The param- 
eters for such a cross section become 

A s =4at 


A = 4 (a+ b) t—pi 


A.— 0 


4 2 nrb 

- — cos — - - 
nir p 


(n even) 


( n odd) 


(A 15) 


With equations (Al5) the parameter N f shown in equation (Al2) becomes 


’9 9 ’9 co COS 1 

AT ^ 2 7^ 2 l 2 p p I,., 

N i — 2 a 7 2 ^ ' / 2 O ^ B 

Sk s 4 fC s O' n = l,3,5 ks 7} 2 1-^2 2 

i "I 


8k s 2 8 k s 2 a 




-+ s 


jZ—i ]„ 2 I 

71 = 1 , 3,5 T) 2 | no2 « = 1 , 3,5 


V BS+n* 


Each of the infinite s umm ations in equation (A17) can now be written in closed form as shown in reference 5, and the 
closed expression for N t then becomes 


1 8k s 2 


K p 

4k s B t d 


cosh 


+ tanhA%i? i K 1 w 
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